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Aims. We aim to understand the star formation associated with the luminous young stellar object (YSO) IRAS 18345-0641 and to 
address the complications arising from unresolved multiplicity in interpreting the observations of massive star-forming regions. 
Methods. New infrared imaging data at sub-arcsec spatial resolution are obtained for IRAS 18345-0641. The new data are used along 
with mid- and far-IR imaging data, and CO (/ = 3-2) spectral line maps downloaded from archives to identify the YSO and study the 
properties of the outflow. Available radiative-transfer models are used to analyze the spectral energy distribution (SED) of the YSO. 
Results. Previous tentative detection of an outflow in the H2 (1-0) SI line (2.122yL/m) is confirmed through new and deeper obser- 
vations. The outflow appears to be associated with a YSO discovered at infrared wavelengths. At high angular resolution, we see 
that the YSO is probably a binary. The CO (3-2) lines also reveal a well defined outflow. Nevertheless, the direction of the outflow 
deduced from the H2 image does not agree with that mapped in CO. In addition, the age of the YSO obtained from the SED analysis 
is far lower than the dynamical time of the outflow. We conclude that this is probably caused by the contributions from a companion. 
High-angular-resolution observations at mid-IR through mm wavelengths are required to properly understand the complex picture of 
the star formation happening in this system, and generally in massive star forming regions, which are located at large distances from 
us. 
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1. Introduction 

It is now understood through recent studies that massive stars, 
at least up to late-0 spectral types, form primarily through disk 
accretion and driving outflows (e.g. Arce et al. 120071 Beuther et 
al.p002b, Davis et al. 2004, Varricatt et al.'20T0^. However, they 
are known to form in clusters in giant molecular clouds that are 
located mostly along the galactic plane at large distances from 
us. The spatial resolution and wavelength coverage of most of 
the studies available are insuflftcient to probe their multiplicity, 
which often lead to inaccuracies in the estimates of the properties 
of the YSOs, their disks and outflows. Detailed studies of indi- 
vidual sources are therefore necessary, which were often not fea- 
sible in surveys though which they were identified. In this paper, 
we present the results of a multi- wavelength observational study 
ofthe luminous YSO IRAS 18345-0641 (hereafter IRAS 18345) 
and discuss the complications associated with treating a YSO as 
single, where it is possibly composed of two or more compo- 
nents. 

IRAS 18345 was identified as a candidate ultra-compact 
Hii region by van der Walt, Gay lard & MacLeod J1995I) , 
who detected a 6.7 GHz methanol maser from its vicinity. 
Highly variable and multi-component 6.7 GHz Class-I methanol 
maser emission has been later detected towards this source by 
several other investigators (Szymczak, Hrynek & Kus i20001 
Szymc zak et al. |20DU1 Sridharan et al. l20Q2l Beuther et al. 
I2002c[ Pestalozzi, Minier & Booth '2005', Bartkiewicz et al. 
120091) . Bronfman, Nyman & May (1996) detected CS(2-1) 
emission from a dense core associated with IRAS 18345 at 
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Vl5/?=95.9 km s"^. A distance of 9.5 kpc has been estimated 
with the distance ambiguity solved (Sridharan et al. 120021) . 

H2O maser emission is used as a tracer of the early stages 
of massive star formation, and is believed to be excited by jets 
(e.g. Felli, Palagi & Tofani 1992). This source hosts H2O maser 
emission (Sridharan et al. 2002; Beuther et al. I2002cl) , showing 
its youth. OH maser emission (1612 MHz), which is known to 
be a tracer of circumstellar disks, has also been detected towards 
IRAS 18345 (Edris, Fuller & Cohen 2007; V^,«;^=93.52 km s'^). 

Beuther et al. (I2002bl) mapped a massive molecular out- 
flow from this source in the ^^CO (/ = 2-1) line. The out- 
flow was seen to be oriented NW-SE, with a projected length 
of 40^X1.84 pc) and a collimation factor /e=1.5. C^^O (2-1) ob- 
servations by Thomas & Fuller (2008 ) detected a line that was 
best fit by a narrow and a broad Gaussian component at 95.65 
and 95.79 km s"^ respectively indicating that they are likely 
to be associated with the same core or material; they suggest 
that the broad component is related to the outflow. Beuther et al. 
(2002a') detected a possible outflow source at 1.2 mm (integrated 
flux density =1.4 Jy) using MAMBO array and the IRAM 30-m 
telescope. This source was also detected at 1 . 1 mm with an in- 
tegrated flux of 0.92(±0.12) Jy in the Bolocam Galactic Plane 
Survey conducted using the Caltech Submillimeter Observatory 
(Rosolowsky et al. 2010), and at 450 and 850yum in the survey 
of high-mass protostellar candidates using JCMT and SCUBA 
(Williams, Fuller & Sridharan [2004]). 

Varricatt et al. (2010) imaged this region in the near-IR K 
band and in a narrow-band filter centred at the wavelength of the 
H2 (1-0) SI line at 2.122yum as a part of their imaging survey of 
massive YSOs. They had a tentative detection of an outflow as an 
H2 line emission knot (MHO 2211; see Figure A20 of Varricatt 
et al). They detected two extremely red near-IR sources (labelled 
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Table 1. Log of WFCAM observations 



UTDate Filters 

(yyyymmdd) used 



Exp. Total int. FWHM 

time (s) time (s) (arcsec) 



20120401 
20120407 
20120426 
20120501 
20120522 



K 

IH 

J 

H2,K,K 



5,1 

10,5 

10 

40 

40,5,1 



180,72 
180,360 
360 
800,800 



0.82,0.84 
1.11,0.99 
0.66 
0.95,0.83 



800,180,72 1.26,1.18,1.16 



'A' [a = 18:37:17.00, S = -06:38:24.5] and 'B' [a = 18:37:16.91, 
S = -06:38:30.7|J£]). Source 'B' detected by them was very close 
to the IRAS source, the 1.2-mm continuum peak of Beuther et 
al. (I2002ab and the CH3OH masers detected by Beuther et al. 
(12002a . 

IRAS 18345 is likely to be in a pre-UCHii phase. Radio 
emission is very weak towards this source. Hofner et al. (2011) 
detected a faint (210yuJy) unresolved 1.3-cm radio contin- 
uum source located at ck= 18:37: 16.91 ^=-06:38:30.4. They pro- 
pose that the radio emission was from an ionized jet driving 
the massive CO outflow detected by Beuther et al. (l20Q2bl) . 
(Sridharan et al. (2002) reported a 27-mJy 3.6-cm source to- 
wards IRAS 18345. However, Hofner et al. (|20TT1) pointed out 
that this emission was from an unrelated source located 3'SE of 
the IRAS source). 



2. Observations and data reduction 

2.1. UKIRTdata 

2.1.1. WFCAM near-IR data 

We observed this region on multiple epochs using the United 
Kingdom Infrared Telescope (UKIRT), Hawaii, and the UKIRT 
Wide Field Camera (WFCAM; Casali et al. l2QQ7l) during the 
UKIDSS backup time. Observations were performed using the 
near-IR /, H and K MKO consortium filters and a narrow 
band filter centered at the wavelength of the H2 (1-0) SI line 
at 2.1218 Jim. WFCAM employs four 2048x2048 HgCdTe 
HawaiiII-2RG arrays; for the observations presented here, we 
used the data from only one of the four arrays in which the ob- 
ject was located. At a pixel scale of 0.4'' pixel" \ each array has 
a field of view of 13.5'xl3.5'. Observations were carried out by 
dithering the object to 9 positions separated by a few arcseconds 
in /, H and K bands, and to 5 positions in H2, and using a 2x2 
microstep, resulting in a pixel scale of 0.2'' pixel" ^ Table]l]gives 
a log of the WFCAM observations performed. 

The data reduction, and archival and distribution were car- 
ried out by the Cambridge Astronomical Survey Unit (CASU) 
and the Wide Field Astronomy Unit (WFAU) respectively. Since 
the photometric calibrations are performed using point sources 
present in all dithered frames and using their 2MASS magni- 
tudes, the photometric quality is good even in the presence of 
clouds. The photometric system and calibration are described in 
Hewett et al. ( 2006 ) and Hodgkin et al. ( 2009 ) respectively. The 
pipeline processing and science archive are described in Irwin 
et al. (2004) and Hambly et al (2008). The left panel of Figure[T] 
shows a JHK colour composite image created from the WFCAM 
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Table 2. Near-IR magnitudes of sources 'A and 'B'. A blank 
space is given when the observations were not done in a particu- 
lar filter on any given UTDate. 



UTdate 
YYYYmmdd 



Magnitudes^ 
H 



K 



Source 'A'; 

20070826^ 

20120407 

20120426 

20120522 



= 18:37:17.016, (5=-6:38:24.36 
not resolved 12.82(0.01) 
not detected 12.81(0.01) 
20.14(0.15) 



Saturated 



8.82 (0.01) 



Source 'B'; Qf=18:37:16.896, (5=-6:38:30.48 



2003 1009^^ 

20070826^ 

20120401 

20120407 

20120522 



not detected 
not detected 



15.934 (0.016) 



16.17(0.05) 



12.69 (0.01) 
12.45 (0.01) 
12.82 (0.01) 

12.53 (0.01) 



18:37:7.00 and 



18:37:6.91 for 'A' and 'B' respectively. 



^The magnitudes are in the UKIRT photometric system, measured us- 
ing an aperture of diameter 2'' and applying aperture corrections. 
The values given in parenthesis are the 1-cr internal errors. ^From 
UK IDSS GPS DR6+. '^From the images presented in Varricatt et 
al. (l20T0l observed using UKIRT and UFTI; the data of all other 
nights were obtained using UKIRT and WFCAM. 



data in a 1.5'xl.5' field centred on IRAS 18345. Sources 'A' and 
'B' are labelled on the figure. Both sources are detected in K. 
This region was also observed in the JHK bands using WFCAM 
as a part of the UKIRT Infrared Deep Sky Survey (UKIDSS; 
Lawrence et al. 2007). We downloaded the JHK images (ob- 
tained for the Galactic Plane Survey, GPS) and merged catalogs 
from the UKIDSS GPS DR6+. 

Source 'A' is saturated in our images with 5 sec exposure, 
so its K magnitude is derived from a mosaic obtained with per 
frame exposure time of 1 sec on 20120522 UT. It is detected well 
in H band. We detect a very faint source at the location of 'A' in 
the / band image obtained on 20120426 UT, when the seeing 
was good. It has a faint neighbour ~0.8"NE. Hence its magni- 
tudes are measured using an aperture of radius 0.4" and applying 
aperture corrections. 'B' is detected well in ^ as a point source 
embedded in nebulosity. Its estimated magnitudes are aff'ected 
by the presence of nebulosity, with the ^-magnitudes in a 3"- 
diameter aperture brighter than that in a 2" -diameter aperture by 
0.21 mag. 'B' is not detected in / and only the nebulosity asso- 
ciated with it is detected in H. Table [2] shows the magnitudes of 
the two sources measured on diff'erent epochs. 

The H2 images were continuum subtracted using scaled K 
band images. The background emission was fitted and removed 
from the H2 and K band mosaics. The seeing obtained in our 
K band mosaics was better than in H2. For any H2 mosaic, the 
K band mosaic with the closest seeing was used for continuum 
subtraction, after smoothing the K band mosaic with a suitable 
Gaussian to match FWHM in the H2 mosaic. The average of 
the ratio of counts ^/H2 was evaluated for a few isolated point 
sources present in both K and H2 mosaics. The background- 
subtracted K mosaics were scaled by these ratios and were sub- 
tracted from the background- subtracted H2 mosaics to perform 
the continuum subtraction. The three H2 mosaics were contin- 
uum subtracted, smoothed with a Gaussian of FWHM=2 pix- 
els to enhance the appearance of the faint line emission features 
and averaged. The right panel of Figure [T] shows the continuum- 
subtracted H2 image of the same region as in the JHK colour 
composite in the left. 
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Fig. 1. The left panel shows WFCAM JHK colour-composite image (/-blue, //^-green, ^-red), and the right panel shows the 
continuum-subtracted H2 image, smoothed with a Gaussian of FWHM=2 pix, and averaged over all our observations. The im- 
ages are constructed from our new WFCAM data, and cover a filed of view of 1.5'xl.5'. The dashed yellow line on the JHK image 
shows the location and direction of the major axis (length x50) of the elliptical distribution of the 6.7 GHz methanol maser spots 
detected by Bartkiewicz et al. (120091 ). The blue (continuous) and red (dashed) contours overlaid on the H2 image show the decon- 
volved blue- and red-shifted lobes of the outflow detected in CO (J=3-2) (see |2.3.1| ); the blue contours are at 3, 5, 10, 15, 20, 25, 
30 and 40 cr, and the red contours are at 3, 5, 7, 9, and 1 1 cr above the mean background. The H2 line emission features are enclosed 
in ellipse and labelled "MHO 2211". The white circle on the top-right corner shows the JCMT beam size (14.6''at 345 GHz). 



2.1 .2. IRCAM3 imaging in L and M' bands 

L and M' imaging observations of IRAS 18345 were per- 
formed using UKIRT and the facility 1-5-yum imager IRCAM^ 
IRCAM3 uses a 256 x 256 InSb array and gives a total field 
of view of 20.8''x20.8''with an image scale of 0.081'' pix'^ at 
the Cassegrain focal plane of UKIRT. The V and M' observa- 
tions were performed on 20000715 and 20000716 UT respec- 
tively. Individual exposure times of 0.2 sec and 0.12 sec with 
100 coadds were used in V and M' respectively. Observations 
were performed by off'setting the field to 4 positions on the array, 
separated by 15''each in RA and Dec. This 4-point jitter pattern 
was repeated thrice. The adjacent frames were subtracted from 
each other and are flat fielded. This resulted in a mosaic with 
two positive and two negative beams for the target. Two of the 
four off'set positions were outside the array and in one of the re- 
maining two, source 'A' was outside the array; hence, the total 
exposure time for 'A' and 'B' were 60 and 120 sec respectively 
in L and 36 and 72 sec respectively in M\ 

The sky conditions were photometric on both nights. The 
excellent seeing on 20000716 (-0.34'') during the M' observa- 
tions enabled us to resolve source 'B ' into two components 'B 1 ' 
and 'B2' for the first time. These two components are at a sep- 
aration of ~0.45"in the SE-NW direction, with 'Bl' located SE 
of 'B2'. Figure |2] shows the M' image of the central region of 
IRAS 18345 with the sources 'A' and 'B' detected. The contours 
generated from the ^-band image are overlaid on Figure |2] The 
insets on the right show a 1.5" x 1.5"field containing source 'B' 
in M'. 'Bl' and 'B2' are labelled on the figure. 



When compared with that in M', the seeing was poor and 
variable during the L' observations (FWHM~0.61"in the mo- 
saic), hence 'B' is not resolved in the final mosaic in L'. 
However, the seeing was better (FWHM~0.38") for two off'- 
set positions (40-s exposure) of the L' observations; an L' mo- 
saic in a 1. 5" X 1.5" field around 'B' constructed using those two 
frames is shown in the bottom-left inset in Figure [2] Photometric 
standards were observed in L' and M' for flux calibration. The 
magnitudes were derived using an aperture of diameter 4.07" (50 
pixels) and applying extinction corrections for the difl'erences in 
the airmass between the object and the standard stars. The M'- 
band photometry of 'B 1 ' and 'B2' were derived using a 0.326" (4 
pixels)-diameter aperture and applying aperture corrections (de- 
rived from source 'A' in the M'-band mosaic, and from standard 
stars for the L' frames). Source 'B' detected in our ^-band image 
is centred on 'Bl', with a nebulous extension in the direction of 
'B2'. Astrometric calibration of the L' and M' images were per- 
formed using the coordinates of 'A' and 'Bl' measured from the 
^-band image. The coordinates and magnitudes of the sources 
detected are given in Table |3] 

2.1 .3. Michelle imaging at 1 2.5yum 

We observed IRAS 18345 using using UKIRT and Michell^on 
20040329 UT. Michelle (Glasse et al. 1997) is a mid-infrared 
imager/spectrometer using an SBRC Si: As 320x240-pixel array 
and operating in the 8-25 yum wavelength regime. It has a field 



^ http ://w w w.j ach . hawaii . edu/UKIRT/instruments/ircam/ircam3 . html 



^http://www.jach.hawaii.edu/UKIRT/instruments/michelle/ 
michelle.html 
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Table 3. IRCAM3 imaging in U and M' bands and Michelle imaging at 12.5 yum 



Source RA"" 

identification 



Dec^ IRCAM3 magnitudes^ 



Michelle flux (Jy)^ 
12.5yL/m 



A 


18:37:17.005 


-06:38:24.5 


5.15(0.01) 


4.59 (0.04) 


4.53(0.18) 


B1+B2 






8.53 (0.02) 


6.71 (0.05) 


1.56(0.10) 


Bl 


18:37:16.917 


-06:38:30.77 


9.65 (0.09) 


7.21 (0.06) 




B2 


18:37:16.897 


-06:38:30.43 


10.12(0.09) 


7.98 (0.06) 





^From the IRCAM images. ^The values given in parenthesis are the 1-cr errors - in the zero points for IRCAM3, and in the flux measured in 
the 4 beams for MICHELLE. 




o 

Q 



17.4 17.2 18:37:17.0 16.8 16.6 

Right asoension 

Fig. 2. M'-band image of the central 15''xl5''field around 
IRAS 18345 Contours from the UFTI ^-band image of Varricatt 
et al. (2010) are overlaid. The bottom-left and bottom-right in- 
sets show 1.5''xl.5''fields surrounding source 'B' in L and M' 
respectively. 



of view of 67.2'' x50.4'' with an image scale of 0.2r7pix at the 
Cassegrain focal plain of UKIRT. 

The Michelle observations were obtained using a filter cen- 
tered at 12.5 yum with 9% passband. The sky conditions were 
photometric. BS 6705 and BS 7525 were used as standard stars. 
The observations were performed by chopping the secondary 
in the N-S direction and nodding the telescope in the E-W by 
15''each (peak-peak) to detect faint sources in the the presence 
of strong background emission. Data reduction was performed 
using the UKIRT pipeline orac-dr (Cavanagh et al. 2008 ) and 
Starlink kappa (Currie et al. 12008b . The resulting mosaic has 
four images of the sources detected (2 positive and 2 negative 
beams). The final image was constructed after negating the neg- 
ative beams and combining all four beams. Figure |3] shows a 
15''xl5''field of the mosaic surrounding IRAS 18345. An aver- 
age FWHM of 0.9 8'' was measured from the sources detected in 
the 12.5 yum image. The astrometric corrections were applied by 
adopting the position of source 'A' from Varricatt et al. (2010). 
The coordinates derived for the fainter 12.5-yum source agrees 
well with that of the near-IR source 'B ' detected by Varricatt et 




17.4 17.2 10:37:17.0 1&.Q 1&.& 

Right asoension 

Fig. 3. Michelle 12.5-yum image of IRAS 18345 covering a field 
of 15''xl5'^ 



al. (12010b . The fluxes measured for 'A' and 'B' are given in Table 
[3] At this spatial resolution, the two components of 'B ' are not 
resolved. 



2.2. Methanol observations 

The J=5-4 (241.80 GHz) and J=7-6 (338.53 GHz) transitions 
of methanol were observed using the 15-m James Clerk Maxwell 
Telescope (JCMT), Hawaii, in July and August 2004. The re- 
ceiver RxA (211-279 GHz) was used for the / = 5 - 4 transi- 
tion and receiver RxB3 (315-373 GHz) was used for J=l-6. In 
both cases the back end was the DAS (Digital Autocorrelator 
Spectrometer) with a bandwidth and spectral resolution of 
250 MHz and 0.1 6 MHz respectively for /=5-4, and 500 MHz 
and 0.625 MHz respectively for /=7-6. Pointing at the JCMT 
was checked every hour and had an rms <2.5''. The weather 
was good during the observations, with opacities in the range 
0.1 3< T225 <0.16. The data were reduced using the data reduc- 
tion package specx (Prestage et al. 120001) . Figure |4] shows the 
methanol spectra. 

CH3OH, is an organic molecule that originates primarily on 
the surface of grains, usually seen towards hot cores. Organic 
molecules such as these appear in the gas phase once suflftcient 
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Fig. 4. CH3OH (/=7-6; 338.53 GHz) and (/=5-4; 241.80 GHz) 
spectra observed using JCMT. 

Table 4. Methanol lines detected 



Line ID 



Rest Freq. 
(GHz) 



(K) 



Radial Width 
Vel. (km/s) 



Amp. /T* 
T;(K) (Kkm/s) 



CH3OH (7-6) 












7i-6i E- 338.344605 


69.61 


95.4 


7.25 


0.16 


1.21 


7o-6oA+ 338.408718 


64.98 


95.4 


5.33 


0.28 


1.60 


72-62 E+ 338.721694 


86.31 


95.5 


9.22 


0.06 


0.61 


CH3OH (5-4) 












5o 4oE+ 241.700168 


46.99 


95.2 


8.03 


0.13 


0.97 


5i-4i E- 241.767247 


39.45 


95.2 


8.88 


0.37 


3.54 


5o-4oA+ 241.791367 


34.82 


95.4 


7.16 


0.53 


4.35 


52-42 E- 241.904643 


59.78 


95.2 


5.11 


0.11 


0.52 



ices have been evaporated from the grain surfaces. This implies 
gas temperatures in excess of 100 K (Rodgers & Charnley '2003 V 
These conditions are achieved either through the localized heat- 
ing of the region immediately surrounding the young star, or in 
the shocked regions produced by a powerful outflow. 

Table |4] shows the rest frequencies, observed radial veloci- 
ties, intensities and line widths of the prominent methanol lines 
detected in our spectra. The line identifications are adopted from 
Anderson, De Lucia & Herbst ([T99Q1) and Maret et al. (120051) . 
The lines are mostly centred at the l'lsr of the core, and have 
large width, suggesting that they originate in a warm region sur- 
rounding the core. 

2.3. Archival data 

IRAS 18345 was detected well in the sky surveys conducted us- 
ing the Wide-field Infrared Survey Explorer (WISE; Wright et 
al. [2QTQj and the AKARI satellite (Murakami et aL,,2007). The 
WISE mission surveyed the sky in four bands - Wl (3.4 yum), 
W2 (4.6yum), W3 (Ufim) and W4 (22jim). At the 6.r-l2.0'' 
spatial resolution of WISE, sources 'A' and 'B' are not resolved; 
however, the excellent positional accuracy (0.15'') of the WISE 
data enables the identification of the source. The location of the 



WISE source is in between 'A' and 'B'; at the shorter wave- 
lengths the centroid is closer to 'A', while at 22yum, it is located 
very close to 'B' (see Figure[5]) confirming that 'B' is the leading 
contributor at longer wavelengths. At the FWHM of -5.7 ''at 9 
and ISyum (Ishih ara et a l.,2010), and 30-40" at 65, 90, 140 and 
160yum (Doi et al. l2Q09l) of the AKARI data, sources 'A' and 'B' 
are not resolved. 

This region was observed in the Spitzer GLIMPSE II sur- 
vey using the Infrared Array Camera (IRAC; Fazio et al. 120041) 
in bands 1-4, centered at 3.6, 4.5, 5.8 and 8.0yum respectively. 
Spitzer also observed this region using the Multiband Imaging 
Photometer for Spitzer (MIPS; Rieke et al. 2004) at 24yum and 
7 1 yum. The reduced Spitzer images and point source photom- 
etry were downloaded from the NASA/IPAC Infrared Science 
Archive. Figure |5] shows a colour composite image of a I'xl' re- 
gion constructed from the Spitzer images at 3.6, 4.5 and 5.8 //m. 
The Spitzer catalog gives magnitudes for source 'A' in Bands 
1-3 only and for source 'B' in Bands 1, 3 and 4 only. We 
performed aperture photometry on the Spitzer images using a 
3.6"-diameter aperture in Bands 1-3 and using a 4.8"-diameter 
aperture in Band 4. Average zero points were derived using a 
few isolated point sources in the field. We derive magnitudes 
5.44+0.05, 4.23+0.05, 3.71+0.05 and 3.35+0.04 for source 'A', 
and 9.13+0.05, 6.68+0.05, 5.85+0.05 and 4.61+0.04 for source 
'B' in bands 1, 2, 3 and 4 respectively. The errors given are 1- 
(T in the estimate of the zero points. IRAS 18345 appears to be 
close to saturation in IRAC bands 1 and 3. Hence we do not use 
the magnitudes derived in these bands in the SED analysis. 

IRAS 18345 is saturated in the Spitzer-MWS 24-jim image, 
but detected well at lljim. At an FWHM of 21 "measured on 
source, 'A' and 'B' are not resolved at 71yum. Using MOPEX, 
we estimate a flux of 71.8 Jy for 'A' and 'B' together at 70yum. 
The locations of the Spitzer source detected at 24jim and lOjum 
are also plotted on Figure [5] 

2.3.1. GO(/=3-2)data 

CO observations were performed using the JCMT in 2008. 
The i2CO(3-2) (345.796 GHz) maps were observed on 
20080421 UT, while the i^CO(3-2) (330.588 GHz) and 
C^^O (3-2) (329.331 GHz) maps were simultaneously observed 
on 20080707 UT. All observations were performed in position- 
switched raster- scan mode and using half array spacing with 
HARP (Buckle et al. 2009) as the front end. The ACSIS auto- 
correlator was used as the back end. ACSIS was setup in dual 
sub-band mode for the ^^CO/C^^O observations giving a band- 
width of 250 MHz and a resolution of 0.06 MHz (0.055 km s"^) 
for each line. The ^^CO data were taken in single sub-band mode 
with a bandwidth of 1000 MHz and a resolution of 0.488 MHz 
(0.42 km s"^). The atmospheric opacity at 225 GHz, measured 
with the CSO dipper, was 0.26 for ^^CO and 0.07-0.08 for 
^^CO/C^^O. The telescope pointing was checked and corrected 
before the observations. The pointing accuracy on source is ex- 
pected to be better than 2". 

The CO data were taken from the JCMT Science Archive 
where they had been reduced using the orac-dr pipeline. The 
pipeline first performs a quality- assurance check on the time- 
series data before trimming and de-spiking. It then applies an 
iterative baseline removal routine before creating the final group 
files. 

Figure |6] shows a 3.3'x3.3'field of the ^^CO spectral cube, 
which shows well defined outflow. Figure|7] shows the same field 
in i^CO and C^^O lines. ^^CO(3-2) is also detected in outflow 
whereas in C^^O, we mostly detect the emission from the central 
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Fig. 5. Spitzer-IRAC colour composite image (blue-band 1, 
green-band2, red-band-3) a rxTfield centred on IRAS 18345. 
The symbols o, >h, X, -h, n and A show the locations of the source 
detected by WISE at 3.4yu m and 2 2yL/m, Spitzer-MIFS at 24yum 
and 70//m, Beuther et al. (l2QQ2aft at 1.2 mm using IRAM 30- 
m telescope amd MAMBO, and and Williams et al. (l2Q04i) at 
850 yum using JCMT and SCUBA, respectively. 



core. Figure fs] shows plots of the integrated spectra of ^^CO, 
i^CO and C^^ in a rxrfield around IRAS 18345. The dotted- 
dashed vertical line shows the rest velocity of the cloud. The 
emission at -110 km s"^ (Figures |6] [t] and [s]) is probably not 
related to the outflow. 



3. Results and discussion 

3.1. The outflows 

Through deeper observations in H2 and K, we confirm the ten- 
tative detection of an H2 outflow by Varricatt et al. (2010). The 
continuum subtracted H2 image in Figure [T] shows the emission 
features labelled MHO 2211, which is composed of two knots 
directed away from 'B'. We derive an outflow angle of 95° east 
of north if these are produced by an outflow from a source at the 
location of 'B ' . 

High angular resolution VLBI observations of the 6.7 GHz 
Class-II methanol maser by Bartkiewicz et al. (2009) revealed 
30 maser spots distributed in the form of an ellipse with semi- 
axes 103 mas and 70 mas respectively. They conclude that the 
6.7 GHz methanol masers form in a disk or torus around the 
massive YSOs, at the interface between the disk and the outflow, 
and that these masers form in a stage before the formation of 
the UCHii region. They estimate a radius of 980 AU for the ring 
of distribution of the maser spots, for a distance of 9.5 kpc. The 
coordinates of the brightest maser spot given by Bartkiewicz et 
al. (^=18:37:16.92106 ^=-06:38:30.5017) agrees with those of 
the infrared source 'B2' within 0.25''. It is also noteworthy that 
the major axis of the ring-like distribution of the maser spots is 
at an angle of 90°, which is nearly perpendicular to the direction 
of orientation of the CO outflow, as it would be expected if they 



form in a disk or torus around the YSO. Assuming a circular 
distribution for the maser spots, we infer an outflow angle of 43° 
with respect to the sky plane. 

The CO (3-2) mapping using JCMT (^ 2.3.1| ) reveal a well- 
defined outflow (Figures [6j [8]). Contours generated from the 
blue- and red- shifted lobes of the outflow are overlaid on the 
H2 image in Figure [T] The contours start at 3 cr above the aver- 
age background level in the integrated blue and red lobes of the 
outflow. The outflow is oriented approximately N-S; we measure 
an angle of 161° E of N for the blue lobe. The outflow appears 
to be centred on source 'B' and is in good agreement with the 
outflow from this source detected in ^^CO (2-1) by Beuther et 
al. (2002b). Figure [T] shows that the direction of this outflow de- 
tected in CO is very difl'erent when compared with an angle of 
95° derived for the H2 outflow. The outflow appears to be at a 
high inclination with respect to the sky plane, as inferred from 
the overlapping blue- and red- shifted lobes (Figure [T]) and the 
high velocity range (Figure [8]); we adopt the angle of inclination 
of 43° with respect to the sky plane inferred from the methanol 
masers for the calculations in this paper. 

The CO spectral line data are used to determine the out- 
flow properties. The velocity-integrated blue- and red- shifted 
spectral images were deconvolved using a 2D Gaussian of 
FWHM=14.58'^ which is the beam size of the JCMT at this 
frequency. At 3-cr level on the blue- and red-shifted images, 
the projected outflow length is 28'', and the width for the blue- 
shifted lobe is 47''. After correcting the observed length for an 
angle of inclination of 43°, we derive a collimation factor of 0.81 
for the outflow. Using a distance of 9.5 kpc, and an average of the 
maximum observed velocities of the blue and red-shifted lobes 
of the outflow of 18 km s"^ we derive the outflow length and 
dynamical time as 1.76 pc and 9.58x10^ years respectively. 

The optical depth in ^^CO (tu) can be obtained from a ratio 
of the antennae temperatures (Tp obtained in ^^CO and ^^CO 
assuming ^2CO/^^CO=89, and using Eq.[l](Garden et al. [T99T1) . 



TlC^CO) l-e- 



1 



T^e^co) 1-^-^13 i_^-^ 



(1) 



Since ^^CO sufl'ers from self absorption near the central ve- 
locity of the line (Figure [5]), we estimated r separately in the 
blue and red wings adopting velocity ranges where both ^^CO 
and ^^CO were detected. We derive values of 18.85 and 15.65 
for T12 in the blue and red wings respectively. 

The ^^CO spectrum exhibits a strong emission component at 
-110 km s"^ (Figures [6] and [s]), which is probably not related 
to the outflow, and is arising from the ambient medium. This 
feature was interpolated out from the spectrum before estimating 
the emission in ^^CO from the blue- and red-shifted lobes of 
the outflow. We used velocity ranges 77-90 km s"^ and 102- 
113 km s"^ respectively for integrating the line wing emission in 
the blue- and red-shifted lobes. 

Following Garden et al. ( I1991K the ^^CO column density is 
derived using the equation 



A^: 



3k 



hBJ{J+\) 



(^e 






87r3V(-/+l)[i-e*^] 



/'" 



dv cm 



(2) 



where B is the rotational constant, yu is the electric dipole mo- 
ment of the molecule, / is the lower level of the transition, and 
Tex is the excitation temperature. Substituting the values for B 
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Fig. 6. ^^CO (J=3-2) spectral cubes covering a 3.3'x3.3'field surrounding IRAS 18345 in a velocity range 77.6 km s"^ to 
113.5 km s"^ centred on the CO-line rest frequency. The left and right panels show the spectral cube in two different orienta- 
tions. The central panel shows the same region as in the left panel, but with the faint emission filtered out to show the outflow more 
clearly. The red square shows the plane of the l'lsr of the core. The contours representing the blue- and red-shifted outflow lobes 
projected on to the sky plane are shown in Fig.[T] 




"^'S^^lfe^o-, 



Fig. 7. ^^CO (J=3-2) and C^^O (J=3-2) spectral cubes with the same spatial and spectral coverages as in Figure [6| The left and 
central panels show the ^^CO cube from two difl'erent orientations and the right panel shows the C^^O cube. The red square shows 
the plane of the l'lsr of the core. 
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Fig. 8. The ^^CO, ^^CO and C^^O (J=3-2) spectra integrated in a 
rxr-box around IRAS 18345. The dotted-dashed vertical line 
at 95.5 km s"^ shows the l'lsr of the core. The emission at 
110 km s"^ is from the ambient medium. 



and /d in Eq. 2, 



U + I; [1 -e^^] -J 

where z; is in km s"^ 

Allowing for beam filling factor and ignoring the eff'ect of 
the cosmic background radiation, 



A^: 



hBJ(J+l) 

2.39x10^4^^^^^ 



(rex + 0.92) 



J rib ' 



Ty dV 



cm 



16.6 (/+1) ^-i^ J r]t(l-e--^) 

. (4) 

Adopting 30 K for Tex (Beuther et al. "2002b) and 10^ for 
H2/CO, Eq. 4 yields 1.37x10^1 cm'^ and 1.03x10^1 cm'^ re- 
spectively for the column density of H2 in the blue- and red- 
shifted lobes. 

The mass in the outflow lobes can be obtained using the re- 
lation m X nH2 X S , where m is the mean atomic weight, which 
is 1.36xmass of the hydrogen molecule, and S is the surface 
area of the outflow lobe (Garden et al. 119911 ). At 3cr above the 
mean background, we measure projected sizes (major axis x mi- 
nor axis) of 55''x45''and 30''x22''respectively for the blue- and 
red-shifted lobes of the outflow. Integrating within these regions, 
within the velocity ranges of 77-90 km s"^ and 102-1 13 km s"^ 
respectively, we get 123 M© and 25 M© for the masses con- 
tained in the blue- and red-shifted lobes of the outflow. This 
gives a total mass of 148 M© for the material in the outflow. 
Using the maximum observed velocities (L'max) 18.5 km s"^ and 
17.5 km s"^ respectively for the blue and red wings, and ignor- 
ing the angle of inclination of the outflow (see Cabrit & Bertout 
1992), we derive an outflow momentum of 2662 M© km s"\ 
energy of 45x10"^^ ergs, and a mass entrainment rate (Mout) of 
15.6x10"^ M© year"^ The total mass in the outflow estimated by 
us is very similar to 143 M© estimated by Beuther et al. (i2002bi) 



from CO (2-1). However, the outflow momentum, energy and 
mass entrainment rate estimated by us are higher. The difl'erence 
comes from the higher value of L'max used by us, and the outflow 
momentum, energy and mass entrainment rate estimated here 
should be treated as upper limits only. 

The collimation factor /c=0.81 we obtain from the C0(3- 
2) map is lower than a value of 1.5 derived by Beuther et al. 
(2002b) from their CO (2-1) map observed at a higher spatial 
resolution. Note that the deconvolved image of the blue- shifted 
lobe is highly asymmetric (see the blue contours in Figure [T]). If 
this asymmetry is due to two overlapping outflows, the collima- 
tion factor will be higher; using a width of 26'' of the red- shifted 
outflow lobe, we get a collimation factor of 1.46. With the di- 
rection of the outflow derived from the aligned H2 knots much 
difl'erent from the direction of the main outflow mapped in CO, 
and roughly agreeing with the direction of the asymmetry of the 
blue lobe of the CO outflow as seen in Figure [T] it is very likely 
that this region hosts at least two YSOs in outflow phase. 

3.2. The sources driving the outflows 

Figure [T] shows that the CO outflow is centred very close to the 
near-IR source 'B'. The H2 line emission knots (MHO 2211) 
also appear to be tracing back to the location of 'B ' . Figure |9] 
shows Si (J - H) - (H - K) colour-colour diagram constructed 
using the objects detected in a 2'x2'field around IRAS 18345. 
Since the UKIDSS observations were performed in better sky 
conditions than our current observations, UKIDSS data are used 
to construct the colour-colour diagram. For 'A' and 'B', the most 
recent magnitudes from Table [2] are used. 'B' for which H and 
/ magnitudes are upper limits are shown with upward and right- 
ward directed arrows to show the directions in which it move in 
the diagram with reliable detection in these bands. Both objects 
are at very large extinction. Source 'B' exhibits large excess and 
is located in the region of the colour-colour diagram occupied by 
luminous YSOs. The ^-magnitudes given in Table [2j observed 
from 2003 to 2012 (all in the same photometric system and using 
the same aperture), show that 'B' is a variable. 

The Spitzer-IRAC colour-colour diagram is useful to identify 
YSOs. Figure [T0| shows the colours of objects detected in the 
IRAC images, within a lO'xlO'field around IRAS 18345. The 
dashed horizontal line shows the boundary between Class I/II 
(below) and Class I objects (above), and the dotted red arrows 
show reddening vectors for Ay = 50 derived using the extinction 
law given in Mathis (1990). Objects with ([3.6]-[4.5])>0.8 and 
([5.8]-[8.0])>l.l are likely to be Class I objects, which are pro- 
tostars with infalling envelopes. Those with 0<([3.6]-[4.5])<0.8 
and 0.4<([5.8]-[8.0])<l.l, within the region shown by the rect- 
angle, are likely to be Class II objects, which are stars with discs. 
Most of the objects are located around ([5.8]-[8.0], [3.6]-[4.5]) = 
0, which are foreground and background stars and discless pre- 
main-sequence (Class III) objects (Allen et al. 2004; Megeath et 
al. [20041 Qiu et al. [20081) . Locations of 'A' and 'B' are shown 
by filled red-circles in Figure [T0| Source "A" is located in the 



region occupied by young sources with discs and large redden- 
ing, and 'B' is located in the region of even younger sources at 
larger reddening and an infalling envelopes. However, note that 
the locations of these objects could be aff'ected by saturation, es- 
pecially in the 8yum band. The locations of the two sources in 
Figures |9] and 10 suggest that 'B' is younger than 'A'. 

All data available with good spatial resolution (UKIRT- 
J,H,K,L\M' and 12.5 yum and Spitzer-IRAC) data show that 
source 'B' has a steeper SED than 'A'. At 12.5 yum, 'B' con- 
tributes 25.6% of the total flux of 'A' and 'B' combined, whereas 
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Fig. 9. WFCAM JHK colour-colour diagram of a 2'x2'region 
surrounding IRAS 18345; UKIDSS data obtained in better sky 
conditions are used. The blue dots show the colours of point 
sources detected. The red dots are the colours of sources 'A' 
and 'B' derived from the most recent observations. The contin- 
uous green lines show the colours of unreddened main sequence 
and red-giant stars taken from Tokunaga et al. J2000I) . The long- 
dashed green line shows the loci of CTTS from Meyer et al. 
(119971) . The long red arrows show the reddening vectors for Ay 
< 50 derived from the interstellar extinction law given by Rieke 
& Lebofsky (1985). 



in U (3.75 yum) and M' (4.7 yum), 'B' contributes only 4.3% and 
12.4% respectively of total flux of 'A' and 'B' combined. At 
shorter wavelengths (UKIRT, Spitzer-IRAC and WISE 3.4yL/m), 
'A' is the leading source. The longer wavelength observations 
(Spitzer-MIFS, WISE-24yum, JCMT and IRAM) do not have 
the spatial resolution to resolve 'A' and 'B'. However, the ex- 
cellent positional accuracies of these observations enable us 
to identify the leading YSO with good confidence. Figure |5] 
shows a colour composite image of a Txl'field surrounding 
IRAS 18345, constructed from Spitzer-IRAC images in bands 1, 
2 and 3. Locations of sources detected in diff'erent observations 
with good positional accuracies are overlaid on the image. As 
can be seen in Figure |5] all longer wavelength observations peak 
near source 'B'. These factors leads us to infer that 'B (Bl-hB2)' 
or some embedded sources very close to it are the YSOs respon- 
sible for the outflows in this system. 

3.2.1 . The spectral energy distribution 

We modelled the SED of IRAS 18345 using the SED fitting tool 
of Robitaille et al. (.2007 j , which uses a grid of 2D radiative 
transfer models presente d in R obitaille et al. (l2QQ6b , developed 
by Whitney et al. (2003 a", '2003b, etc). The grid consists of SEDs 
of 20000 YSO models covering a range of stellar masses from 
0. 1 to 50 M© and and evolutionary stages from the early envelope 
infall stage to the late disc-only stage, each at 10 diff'erent view- 



Fig. 10. The IRAC colour-colour diagram. The open (black) cir- 
cles show the objects detected in all four IRAC bands in a 
lO'xlO' field centred on IRAS 18345. The locations of sources 
'A" and "B" are shown by filled (red) circles. The loci of Class II 
sources are shown by the blue rectangle; the dashed horizontal 
line shows the boundary between Class I/II objects (below) and 
Class I objects (above). The dotted red arrows show reddening 
vectors for Ay=50. 



ing angles. With all the longer wavelength data peaking close to 
'B', we have treated that as the leading YSO in this region. 

Up to 12.5yum, we have fluxes of 'B' available. WFCAM K, 
IRCAM UM\ Spitzer-IRAC Bands 1 and 3 and MICHELLE 
12.5-yum flux of 'B' were used. WFCAM H magnitude was not 
used since in H, we are probably seeing only the nebulosity as- 
sociated with 'B'. Spitzer-IRAC Bands 2 and 4 fluxes were not 
used since they are near saturation. The resolution of WISE im- 
ages is not sufficient to resolve the two sources, so WISE 3.4- 
12-yum data were not used to construct the SED. At 22yum, the 
WISE source is located close to 'B ' , there is still likely to be a 
significant contribution to the flux from 'A'. Hence the 22-jim 
flux was used as an upper limit. IRAS catalog lists the flux val- 
ues for IRAS 18345 in 12 and lOOyum bands as lower limits, 
and the 100-yum flux may be afl'ected by infrared cirrus, so those 
values were not used in the SED analysis. AKARI-IRC 9-yum 
flux was not used since that wavelength was already covered 
by the higher resolution UKIRT data. Similarly, MSX data in 
Bands A, C and D were not used. AKARI-IRC 18-yL/m data and 
MSX Band-E data are treated as upper limits only since they 
contain contributions from source 'A' . The large bandwidths of 
the IRAS and AKARI photometry warrant colour corrections, 
which were performed using the correction factors given in their 
respective point source catalogs (Beichman et al. 1988; Kataza 
et al. 2010; Yamamura et al. GOTOl) . The JCMT-SCUBA 450 and 
850-/im fluxes were adopted from Williams, Fuller & Sridharan 
(120041) . Flux estimates at 1.1 mm and 1.2 mm were adopted 
from Rosolowsky et al. (120101 and Beuther et al. (2002a) re- 
spectively. Figure 



1 1 shows the model fit to the spectral energy 
distribution. Table 5]shows the results of the SED modelling. The 
model fit to the SED yields a very young object of mass 19.8 M©, 
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Table 5. Results from the SED fitting 



10" 



Parameter 


Best-fit values'" 




Stellar mass (M©) 


19.8 (+3, -1.5) 


_ 10-^ 


Stellar age (yr) 


3.3(+.49,-0.77)xl03 


\ 


Stellar radius (R©) 


230 (+16, -20) 


E 


Stellar temperature (K) 


4480 (+260, -167) 




Disk mass (M©) 


1.9(+6,-0)xl0-i 


g.io-^° 


Disk accretion rate (M©yr~^) 


2.5(+4, -1.4)xl0-5 


(U 


Disk/envelope inner radius (AU) 


9.9 (+0,-1.3) 


ur^ 


Disk outer radius (AU) 


11.7 (+11.8,-0) 


r< 


Envelope mass (M©) 


2.2(+3.6,-0)xl03 


10-^^ 


Envelope accretion rate (M©yr~^) 


4.1(+4.3,-0.2)xl0-3 




Envelope outer radius (AU) 


1 (+0, -0.4)xl05 




Total Luminosity (L©) 


1.9(+2, -1.5)xl04 




Ay circumstellar (mag) 


804 (+508, -453) 


10 ^2 


Ay interstellar (mag) 


25 (+13) 




Angle of inclination of the disk axis (°) 


18 





The values given in parenthesis are the rms of the differences, in 
parameter values of models with (x^ ~ xlestfu^ P^^ ^^^^ point < 3, 
above and below those of the best-fit model, estimated with respect 
to the parameter values of the best-fit model. A distance of 9.5 kpc 
is adopted. 



age 3.3x10^ years and luminosity l.QxlO'^Lo. Figure 11 shows 
that there is a large disagreement between the fluxes obtained 
by difl'erent telescopes in the wavelength range of 60-140 yum. 
The flux estimated from the Spitzer-M\PS lOjim image is much 
lower than what was seen in IRAS 60yum and AKARI 66yum 
bands. Similarly, the AKARI 90 and 140-yum fluxes also devi- 
ate significantly from the model SED. Inadequate spatial resolu- 
tion and the difl'erent contributions from the background could 
be severely afl'ecting the fluxes measured in these bands used by 
telescopes of difl'erent sizes and spatial resolutions. The SED 
analysis for this source available in the literature, performed 
using the same model (Tanti, Roy & Duorah, 2011) using all 
the available data, obtained a 12.99 M© YSO with an age of 
1.36x10^ years and a luminosity of 7.21x10^ L©. The difl'erence 
mostly comes from the fact that they do not attempt to separate 
the fluxes of 'A' and 'B'. 

The results obtained from our SED analysis bring some 
complications. The dynamical time scale of the outflow 
(9.58x10^ years) derived by us from the CO (3-2) data is much 
larger than the age of the YSO (3.3x10^ years) obtained from the 
SED analysis. This is unrealistic. This discrepancy may mostly 
be caused by the contribution from a companion at longer wave- 
lengths, thereby making the star appear younger. We have re- 
moved the contributions from source 'A' to the total flux for 
wavelengths up to 12.5 yum. The centroids of all detections at 
longer wavelengths are located close to 'B', hence 'A' is not 
likely to be the main contributor at longer wavelengths. The 
more probable chance is that 'B ' itself is a binary or multiple as 
strongly suggested by its resolution into two components 'B 1 ' 
and 'B2' in U and M\ and the detection of outflows with two 
difl'erent directions from H2 and CO observations, both of which 
appear to be centred near 'B'. High-angular-resolution observa- 
tions at longer wavelengths are therefore required to properly un- 
derstand the complexity of the star formation taking place here. 




1000 



A (//m) 



Fig. 11. The SED of IRAS 18345. The filled circles show the 
data in WECAM K band, IRCAM U and M' bands, Spitzer- 
IRAC Bands 1 and 3, Michelle 12.5 yum, IRAS 60yum, AKARI 
at 65-140 yum, Spitzer-MIFS 70yum, SCUBA 450 and 850yum, 
Bolocam 269 GHz and IRAM 240 GHz data. The downward- 
directed triangles show AKARI-IRC ISyum, MSX band-E and 
WISE band-4 data which are treated as upper limits. The contin- 
uous line shows the best fit model. The grey lines show subse- 
quent good fits for (x^ ~ xlestfit^ P^^ ^^^^ point < 3. The dashed 
line corresponds to the stellar photosphere for the central source 
of the best fitting model, as it would look in the absence of cir- 
cumstellar dust (but including interstellar extinction). 

4. Conclusions 

1. IRAS 18345 is a very young massive YSO in a phase of ac- 
tive accretion. 

2. Through near and mid-IR observations, we identify the cen- 
tral source, which is highly obscured and is probably double 
at a separation of -0.45''. 

3. We confirm the 2.122//m H2 outflow tentatively detected by 
Varricatt et al. (120 101) . The direction of this outflow is very 
difl'erent from the direction of the main outflow mapped in 
CO. The origins of both H2 and CO outflows appear to trace 
back to the location of the infrared source mentioned above. 
In addition, the SED analysis does not fit well with a sin- 
gle YSO model. This suggests that more than one YSO are 
driving outflows here. 

4. The SED analysis yields a 19.8 M© star of age 
3.3 X 10^ years. The age is far lower than what is im- 
plied by the dynamical time of the outflow. The unresolved 
binarity is probably the factor afl'ecting the SED analysis, 
resulting in an overestimate of the mass and underestimate 
of the age. 

5. The collimation factor of the outflow derived by us from the 
CO (3-2) maps is low (0.81). It needs to be explored if this 
is real, or is due to the presence of more than one outflow. 

6. We emphasize the need high-angular-resolution observations 
at mid-IR and longer wavelengths to properly understand 
this system, and in general massive star forming regions. 
High angular resolution interferometric CO line observations 
are also required to understand if the apparently single out- 
flow seen in CO is composed of more than one component. 
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